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First-principles calculations were carried out to investigate interfacial strain effects on spin in- 
jection and spin polarization of a magnetic tunnel junction consisting of half- metallic full-Heusler 
alloy Co2CrAl and ferroelectric perovskite NaNbOs- Spin-dependent coherent tunneling was calcu- 
lated within the framework of non-equilibrium Green's function technique. Both spin polarization 
and tunnel magnetoresistance (TMR) are affected by the interfacial strain but their responses to 
compressive and tensile strains are different. Spin polarization across the interface is fully preserved 
under a compressive strain due to stronger coupling between interfacial atoms, whereas a tensile 
strain significantly enhances interface states and lead to substantial drops in spin polarization and 
TMR. 



I. INTRODUCTION 

Injection of spin-polarized electrons or holes from one 
material to another has been the subject of intensive 
studies because it is essential in many spintronic device 
applications such as magnetic tunneling junction (MTJ), 
spin valves and spin injectors [l|-[3|. In general, there are 
two major approaches to achieve a high spin injection 
efficiency. One is to use a semiconducting or insulating 
barrier to selectively filter carriers with one spin state 
from normal ferromagnetic metals such as Co or Fe, as 
demonstrated in the Fe/MgO/Fe MTjIl. The other is 
to adopt a half metal as the spin source (sf-Q- For the lat- 
ter approach, since there exists a gap for one spin state, 
the spin polarization (SP) at the Fermi energy can reach 
100% and a high spin injection efficiency can be achieved 
without the need of a spin filter. Besides the high spin 
polarization, another significant advantage of using half 
metals instead of ordinary metallic ferromagnets as a spin 
source is the elimination of impedance mismatch between 
the spin source and the semiconductor which is essential 
for efficient spin injection [8^. 

Among the various half metals proposed for spin 
source, the full-Heusler alloys (A2BC) in the L2i struc- 
ture have attracted great interest due to their large mag- 
netic moments per unit celljil, 0- In this type of al- 
loys, there exists a band gap for minority spin arising 
from hybridization between the orbitals of A and 
B atoms and interaction between the d orbitals of the 
nearest neighbor type A atoms Extensive theoreti- 
cal and experimental studies have been carried out to 
explore high SP in several full-Heusler alloys j6!-[9|. For 
example, a spin-injection efficiency of around 50% was 
achieved from a full-Heusler alloy, Co2FeSi, to (Al,Ga)As 
and MgO [lo|, [Hi • A relatively large tunnel magnetoresis- 
tance (TMR) has been reported for Co-based full-Heusler 
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alloys such as Co9Cri _^Fe^ Alj?!, fT2| . 

The measured SP of bulk Heusler alloys, however, is 
usually much lower than theoretical predictions. Most of 
bulk Heusler allo ys l ose their high SP due to phase dis- 
order and defects [l3[. More importantly, when a Heusler 
alloy is heterogeneously grown on another material, inter- 
face states arising from charge-discontinuity at the inter- 
face dramatically reduces its SP and TMR ratio [14]. In 
addition, SP could be affected by interfacial strain due 
to lattice mismatch between the Heusler alloy and the 
semiconductor/insulator. Theoretical calculations have 
shown that strain can alter the width of the minority 
band gap and the position of the Fermi level in the 
minority band gap[15]. Botters et a/.[l6[ demonstrated 
that a small strain (0.08%) applied to half-Heusler alloy 
NiMnSb induces a large change (20%) in its magnetic 
anisotropy. A similar strain-induced effect was also ob- 
served in Co2MnGa[17]. On the other hand, it is much 
harder to control the interfacial strain experimentally due 
to complex film growth/deposition processes. Relatively 
little is known about interfacial strain effect besides it 
varies depending on the thickness of the film[l8[. The 
performance of a spintronic device based on heterostruc- 
tures is dominated by interface quality. While the qual- 
ity of an interface can be due to many factors, interfacial 
strain is critical and often the cause of other interface im- 
perfections such as interfacial defects. It is therefore crit- 
ically important to systematically investigate and quan- 
titatively understand the effect of interfacial strain on SP 
and TMR of Heusler alloy /semiconductor j unctions [l9|- 
121]. 

In the present work, we evaluate effect of an interfacial 
strain on spin injection of a Heusler alloy-based junc- 
tion, CosCrAl/NaNbOa/CosCrAl (CCA/NNO/CCA), 
using a first-principles computational approach. We 
choose Co2CrAl because of its remarkable magnetore- 
sistive property and persistence of its half-metallic be- 
havior at surface and interface 0, |22|. On the other 
hand, sodium niobate is a suitable spacer layer because of 
its broad ferroelectric and piezoelectric applications and 
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good lattice match with Co2CrAl. We demonstrate that 
both SP and TMR are affected by the strain but their 
responses to positive (tensile) and negative (compressive) 
strains are different. 



II. CALCULATION DETAILS 

First-principles calculations based on spin-polarized 
density functional theory were carried out using the 
VASP package [2 3j. The projector augmented wave 
method was used to describe electron-ion interaction 
while the Perdew-Burke-Ernzerhof generalized gradient 
approximations (PAW-PBE) was adopted for electron ex- 
change and correlation. A cutoff energy of 400 eV and a 
6x6x1 Monkhorst-Pack k-mesh were chosen which pro- 
vide a good balance between accuracy and computational 
time. Structural relaxations were performed until the 
Hellmann-Feynman force on each atom became less than 
0.02 eV/A. The CCA(001)/NNO(001) heterostructure 
is formed by directly stacking CCA(OOl) on NNO(OOl) 
without rotation, as shown in Fig. 1. There is a small 
lattice mismatch of 2.4% between CCA and NNO. We 
assume that the Heusler alloy is grown on the top of the 
semiconductor [lo|, and therefore, reduced slightly the in- 
plane lattice constant of cubic CCA to match those of 
NNO, as it is often done in interface modeling [8i]. The 
lattice parameter of the supercell in the direction normal 
to the interface was adjusted until the total energy of the 
system reached a minimum, which resulted in a c/a of 
1.05 for CCA. Odd numbers of CCA (19) and NNO (9) 
atomic planes were used in the model to prevent forma- 
tion of artificial electrostatic fields across the interfaces. 



III. RESULTS AND DISCUSSION 

Since CCA consists of alternate Co and CrAl planes in 
the [001] direction, and similarly NNO is composed of al- 
ternate Nb02 and NaO planes, a CCA(001)/NNO(001) 
interface can be built with Co- or CrAl-terminated CCA 
and Nb02- or NaO-terminated NNO. In addition to the 
four interface models (models 1-4 in Table I) built with 
various combinations of CCA(OOl) and NNO(OOl) ter- 
minations, we consider another model (model 5) which 
is formed by CrAl-terminated CCA and NaO-terminated 
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FIG. 1: (Color online) Atomic structure of the 
CCA/NNO/CCA MTJ. The models are built by stack- 
ing CCA directly on NNO along the [001] crystalline 
direction. 



NNO, but the CCA is shifted related to NNO in the inter- 
face plane so that Cr forms bond with O and Al is bonded 
to Na, instead of the Al-0 and Cr-Na bonds in the in- 
terface built by directly stacking CCA on NNO (model 
4) . To evaluate the relative stability of the various inter- 
face models, we calculate the interfacial cohesive energy 
(Eice) which is defined as 



(£^ccA + Enno - Ecca/nno/cca)/2A (1) 



where Eqca, ^'nno, and £^cca/nno/cca are total 
energies of the CCA slab, the NNO slab, and the 
CCA/NNO/CCA supercell, respectively. 2 A accounts 
for the total area of the 2 interfaces in the supercell. 
The calculated interfacial cohesive energies, along with 
the lengths of atomic bonds at the interfaces, are listed 
in Table L The positive and significant values of the in- 
terfacial cohesive energies indicate that all five interfaces 
are stable which is consistent with the fact that Heusler 
alloys can form stabile junctions with distinct interfacial 
structures. It is interesting to note that both CrAl-NaO 
interfaces are surprisingly more stable than the CrAl- 
Nb02 interface which has double the number of interfa- 
cial oxygen bonds. The strongest CCA/NNO interface 
is the Co-Nb02 interface with O atoms lying on top of 
Co atoms as shown in Fig. 1. The large cohesive energy 
of this interface (2.46 J/m^) is accompanied by a strong 
Co-0 bond. The Co-0 bond length (1.94 A) is close to 
that in bulk CoO. Since the binding energy of this in- 
terface is 20% larger than that of the next most stable 
interface, this model is used in subsequent calculations 
of electronic structure, SP and spin injection. 



TABLE I: Calculated interfacial cohesive energies (Wice) and 
interface bond lengths (R) for the five CCA/NNO interfaces 
being studied. 



Interface 


V^ice (J/m^) 


R{A) 


1: 


Co-Nb02 


2.46 


Co-O (1.94) Co-Nb (2.71) 


2: 


Co-NaO 


1.72 


Co-O (1.85) Co-Na (2.68) 


3: 


CrAl-Nb02 


1.64 


Al-O (1.99) Cr-O (1.97) 


4: 


CrAl-NaO 


2.01 


Al-O (1.83) Cr-Na (2.81) 


5: 


CrAl-NaO 


1.93 


Cr-O (1.89) Al-Na (2.91) 



The spin-polarized local density of states (LDOS) of 
the CCA/NNO/CCA MTJ is shown in Fig. 2. The strain 
introduced by matching the lateral lattice constant of 
CCA to those of NNO induces minor changes in the mag- 
netic property of CCA, compared to that of bulk cubic 
CCA. The calculated local magnetic moment of Co is 
0.73 jiiB and that of orthorhombic Cr is 1.61 /i^, which 
are very close to those in cubic CCA (0.77 /i^ for Co 
and 1.54 fiB for Cr), and the total magnetic moment is 
in good agreement with prediction (3 /i^) based on the 
Slater-Pauling scheme, Mf = — 24, where Mt is the to- 
tal magnetic moment (in /hb) and Zt is the total number 
of valence electrons in the unit cell. 

The most important feature of the LDOS is that the 
half-metalhc property of CCA is preserved at the inter- 
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FIG. 2: (Color online) LDOS for the majority (top panels) 
and minority (bottom panels) spin electrons projected to the 
O, Co and Nb atoms at the Co-Nb02 interface and in bulk 
region, respectively. The vertical dashed line indicates the 
position of the Fermi level. 



face, despite of substantial modifications of the electronic 
structures of the interfacial atoms compared to those in 
the bulk. For instance, the magnetic moments of the in- 
terfacial Co (0.39 fis) and Cr (1.11 fis) become smaller 
than the corresponding values in bulk CCA. The minor- 
ity band gap is reduced from ~0.5 eV (GGA value) in 
bulk CCA to ~0.2 eV due to orbital mixing between 
the Co 3d and O 2p orbitals. This strong hybridization 
pushes down the bonding state of the interfacial O atoms 
by nearly 1 eV (Fig. 2d). In addition to these bonding 
states between O and Co atoms in the low energy part 
of the 3d states of the interfacial Co atoms, the oxy- 
gen LDOS displays a broad band of antibonding states, 
ranging from —3 eV to 3 eV, which are absent in bulk 
NNO. Similar metal induced gap states (MIGS) appear 
more prominently for the interfacial Nb atoms as shown 
in Fig. 2f. Since the length of the Nb-Co bond is 2.71 A 
and the Wigner radii of Nb and Co atoms are 1.27 and 
1.30 A respectively, direct bonding between the interfa- 
cial Nb and Co atoms seems unlikely, and these MIGS on 
Nb atoms must arise from an indirect coupling between 
Nb and Co mediated by O atom. This indirect coupling 
leads to a substantial exchange-splitting of the states in 
the energy range from the Fermi energy down to about 2 



eV below (Fig. 2f), and a magnetic moment of 0.05 /i^ on 
the interfacial Nb atom which is aligned antiparallel to 
those of the Co atoms. This situation is quite similar to 
the case of Co/SrTiOs junction where the interfacial Ti 
atom is found to possess a small magnetic moment which 
is antiferromagnetically coupled to that of the interfacial 
Co atom [23. 
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FIG. 3: (Color online) Variation of the planar averaged SP 
along the interface normal direction of the CCA/NNO/CCA 
MTJ, at various strains. The position of the Co-Nb02 in- 
terface plane is indicated by the vertical dashed line. The 
horizontal dashed line corresponds to 100% SP in each case. 

Next we discuss the strain effects on the SP and elec- 
tronic tunneling of the MTJ. A uniaxial strain along the 
normal direction of the interface is considered and is re- 
alized in our model through changing the length of the 
supercell while fixing the atomic positions in the extreme 
left and extreme right unit cells of CCA in the supercell 
(Fig. 1). Contraction or expansion in the lateral direc- 
tions are estimated to be small in the range of unaxial 
strain in the present study and are neglected. The sys- 
tem is then relaxed under the constraint of fixed supercell 
parameter c. The SP of each atom at the Fermi level is 
defined as 

SP = (2) 

where and are the majority and minority LDOS, 
respectively, at the Fermi level, projected on the con- 
cerned atoms. The calculated in-plane averaged SP is 
shown in Fig. 3 for various strain values. Since the struc- 
ture and SP are nearly symmetric with respect to the 
middle of the MTJ structure, only SP for half of the 
MTJ is shown in the figure. A 100% SP can be ob- 
served for atoms at the CCA side, reflecting the half- 
metallic character of CCA. At zero strain, an extremely 
high SP (nearly 100%) is seen at the Co-Nb02 interface 
which extends to NNO. Under a compressive strain, a 
steady increase of SP in the NNO part can be clearly ob- 
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served. This is due to stronger coupling between interfa- 
cial atoms under a compressive strain. Since the minority 
band gap of CCA originates from the Cu-hu sphtting due 
to interaction between the nearest neighbor Co atoms 
the sphtting is enhanced under a compressive strain due 
to increased hybridization. Under a moderate positive 
strain, SP in NNO drops shghtly while the polarization 
in CCA remains high (>90%). However, when the strain 
approaches ~4%, the SP decreases dramatically. This 
cannot be attributed to structural instabilities since no 
significant structural changes occur in the heterostruc- 
ture at this strain and the Co-0 bond length increases 
only slightly, from 1.94 A at zero strain to 1.95 A at 4% 
strain. 



transmission of the majority states and much smaller and 
spiky transmission peaks for the minority states suggest 
that a high TMR ratio can be expected. The calculated 
TMR ratio shown in Fig. 4c confirms this prediction. 
Furthermore, the TMR ratio fluctuates with compres- 
sive strain, but decreases rapidly with the tensile strain. 
At 4% of tensile strain, the TMR ratio drops by two 
orders of magnitude compared to the strain free value. 
This result is in good agreement with recent experimen- 
tal observation [28] that the TMR of a Co2MnSi-based 
MTJ fluctuates slightly, without any abrupt changes, in 
the process of applying a hydrostatic pressure on the sys- 
tem. To date, no experiment has been conducted to in- 
vestigate the effect of a positive strain on such a system. 
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FIG. 4: (Color online) In-plane wave vector k// = {kx,ky) 
dependence of majority spin (a) and minority spin (b) trans- 
missions at the Fermi level of the CCA/NNO/CCA MTJ in 
the parallel magnetization configuration without strain, (c) 
TMR ratio as a function of strain. 

According to Jullierejlsj, SP is a key quantity in deter- 
mining TMR of a MTJ. The asymmetrical response of SP 
to positive and negative strains in the MTJ prompted us 
to calculate the electron tunneling and TMR across the 
junction [26]. The TMR ratio is defined as 

G^^ — G^i 



TMR: 



a 



n 



where G^^ and G^^ are conductances of the MTJ when 
the magnetizations of the two electrodes are parallel and 
antiparallel, respectively. The conductance per unit cell 
area is given by the Landauer-Biittiker formula 



where T^jiku) is the transmission probability of the elec- 
tron with spin a at the Fermi energy, with /c^ = (kx^ky) 
being the Bloch wave vector in the plane of the junction. 
The calculated /c^-resolved transmission in the parallel 
magnetization configuration is shown in Fig. 4. Simi- 
lar to the case of the CCA/MgO system [27|, no trans- 
mission was found for the majority spin states around 
the F point. For the minority spin states, the transmis- 
sion shows spiky structures, but of much smaller mag- 
nitudes, which can be attributed to the resonant tun- 
neling through the weak interface states. The strong 




FIG. 5: (Color online) The differences between majority spin 
transmissions at the Fermi level in the parallel magnetization 
configuration in strained and strain free CCA/NNO/CCA 
MTJ: (a) ^^(-4%) - G'tt(O), (c) ^^(4%) - 0^(0). The 
transmission at zero strain is given in (b). The Fermi sur- 
faces corresponding to the two folded bands in the Brillouin 
zone of the orthorhombic Co2CrAl are shown in (d) and (e), 
respectively. 

The features of the transmission spectrum are related 
to the shape of the Fermi surface folded into the first Bril- 
louin zone of the deformed orthorhombic CCA supercell. 
The difference between the /c^-resolved majority trans- 
missions of the strained MTJ and the strain free MTJ in 
the parallel magnetic configuration is shown in Fig. 5a for 
—4% strain and in Fig. 5c for 4% strain, respectively. At 
a strain of —4%, large changes are seen at the highly 
conducting /c-points compared to transmission at zero 
strain which is shown in Fig. 5b. These /c-dependent 
changes can be qualitatively understood by variations of 
tunneling probability of electrons at the Fermi surface. In 
contrast, at 4% strain, new transmission channels, which 
are absent in the transmission spectrum of the strain free 
MTJ, are found in the central region of the Brillouin zone 
where a hole appears in the Fermi surface (Figs. 5d and 
5e). These additional conducting channels are related to 
the interface states resulting from hybridization between 
atomic orbitals of interface atoms due to the strain in- 
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duced interface structural changes. The interface states 
also greatly enhance the conductance of the minority 
electrons and destroy the half- metallic character of CCA 
at the interface. Since the conductance in the antiparallel 
magnetic configuration of the electrodes G^i{kf/) can be 
approximated by the geometrical mean of the majority 
and minority conductances obtained for the parallel con- 
figuration, i.e., G^i{k//) = ^/G^^{k//)Gll{k//), one can 
expect G^i{k//) to increase dramatically for states around 
r and a substantial reduction in the TMR at 4% strain, 
as shown in Fig. 4c. 



IV. CONCLUSIONS 

In conclusion, results of our first-principles calculations 
show that both SP and TMR of the CCA/NNO/CCA 



MT J present asymmetrical responses to positive and neg- 
ative interfacial strains. While a high SP is maintained at 
the interface and the TMR ratio fluctuates slightly under 
compressive strain, a positive strain of 4% leads to a drop 
of TMR by ^2 orders of magnitude due to electron tun- 
nehng through the interface states. Our study shows that 
a strain in the MTJ can result in significant changes in 
transport properties of the Heusler /semiconductor junc- 
tion due to strain induced changes in coupling of elec- 
tronic states across the interface. This study provides a 
direct evaluation of the effects of interfacial strain devel- 
oped in the films on spin injection of Heusler alloy-based 
spintronics devices. 
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